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ABSTRACT 

Recent H I surveys of the Large Magellanic Cloud (LMC) with the Australia 
Telescope Compact Array (Kim et al. 1998) and the Parkes multibeam re- 
ceiver (Staveley-Smith et al. 2003) have focussed, respectively, on the small-scale 
(< 20') structure of the interstellar medium (ISM) and the large-scale (> 1°) 
structure of the galaxy. Using a Fourier-plane technique, we have merged both 
data sets providing an accurate set of images of the LMC sensitive to structure 
on scales of 15 pc (for an LMC distance of 50 kpc) upwards. The spatial dynamic 
range (2.8 orders of magnitude), velocity resolution (1.649 km s _1 ), brightness 
temperature sensitivity (2.4 K) and column density sensitivity (8.9 x 10 18 cm -2 
per 1.649 km s _1 channel) allow for the studies of phenomena ranging from the 
galaxy-wide interaction of the LMC with its close neighbours to the small-scale 
injection of energy from supernovae and stellar associations into the ISM of the 
LMC. This paper presents the merged data and the size spectrum of H I clouds 
which is similar to the typical size spectrum of the holes and shells in the Hi 
distribution. 

Subject headings: galaxies: individual (Large Magellanic Cloud) — galaxies: ISM 
- ISM: neutral hydrogen — Magellanic Clouds — radio lines: galaxies 
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1. Introduction 

High resolution H I observations of nearby galaxies allow the study of many aspects 
of their dynamics, morphology and interstellar physics. Measurement of the vertical Hi 
distribution and velocity dispersion, for example, allows the mass-to-light ratio and dark 
matter content of galaxy disks to be studied (Oiling 1996). The shape of the inner velocity 
field, using Hi as a tracer, can discriminate between different halo models (de Blok et 
al. 2001). Also, the relationship between H II regions and Hi holes says much about the 
evolution of galaxies and the propagation of star-forming regions (Walter & Brinks 2001). 

The LMC, as the nearest gas-rich galaxy to our own, has been the subject of several H I 
studies (Luks & Rohlfs 1992; McGee & Milton 1966). The advantage of studying the LMC 
in Hi is that it is the nearest galaxy to our own with a distance of 50—55 kpc (Feast 1991), 
it is presented nearly face-on, and it is a very gas-rich and active star forming galaxy, so that 
it allows for a detailed study of the structure, dynamics and the interstellar medium of a 
star-forming galaxy at close range. The early single-dish observations revealed the existence 
of expanding Hi shells in the LMC (Westerlund & Mathewson 1996; McGee & Milton 1966). 
A very spectacular example of such a shell associated with the major star-forming region 
Constellation III was described by Dopita, Mathewson, & Ford (1985). 

However, its location in the far south means that high resolution studies had to await 
the advent of the Australia Telescope Compact Array 1 (ATCA). The advent of ATCA finally 
allowed us to take full advantage of the aperture synthesis technique in southern hemisphere 
observations, and our recent high-resolution H I survey of the Large Magellanic Cloud (LMC) 
revealed that the structure of the neutral atomic interstellar gas to be dominated by numerous 
shells and holes as well as complex filamentary structure (Kim et al. 1998). Its huge angular 
scale (~8° for the inner disk), required effective observing techniques and mosaicing methods 
(Staveley-Smith et al. 1997), and required the development of new deconvolution algorithms 
(Sault, Staveley-Smith, & Brouw 1996). 

At the largest scales, understanding the interaction of the LMC with the Galaxy and 
the Small Magellanic Cloud (SMC) is important (Putman et al. 1998; Weinberg 2000). 
Also important is knowing the viewing geometry of the disk of the LMC (van der Marel & 
Cioni 2001), its mass and dark matter content (Kim et al. 1998; Alves & Nelson 2000) and 
the origin of the off-center bar (Zhao & Evans 2000; Gardiner, Turfus, & Putman 1998). 
At intermediate scales, the origin of the supergiant shells is tantalizing. Are supernovae 
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and stellar winds (Dopita, Mathewson, & Ford 1985; Chu & Mac Low 1990) sufficient, or 
are other factors such as instabilities in the ISM (Wada et al. 2000), high- velocity cloud 
collisions ( Tenor io-Tagle & Bodenheimer 1980), or gamma ray bursts (Efremov, Ehlerova, 
& Palous 1999) important? At smaller scales, the structure of the ISM, the structure of 
photo-dissociation regions and the detailed feeding and feedback of star-forming regions are 
important for understanding their evolution. 

In this paper, we perform the last step of our high-resolution Hi survey of the LMC 
by combining the already-made ATCA observations (Kim et al. 1998) with the relatively 
newer observations made with the multibeam receiver (Staveley-Smith et al. 2002) on the 
Parkes telescope. This step is necessary in order to provide a single-dish data cube containing 
structure over the full range of spatial frequencies on the sky. The ATCA images, on their 
own, may not provide reliable Hi masses for objects extended over more than 10' - 20' as 
the minimum ATCA baseline for these observations was 30m, for which the fringe spacing 
is X/D = 24' (though mosaicing extends the effective fringe spacing to ~ 36' - see Holdaway 
1999). Similarly, Parkes Hi images have a resolution of 14'- 16', depending on the details 
of the gridding of the data, and are therefore not useful for examining structures smaller 
than ~ 25'. Several examples of combination of aperture synthesis and single-dish data can 
be found in the literature: Ye & Turtle (1991) and Stanimirovic et al. (1999) utilize image- 
plane combination, whereas Bajaja & van Albada (1979), and Wilner & Welch (1994) utilize 
Fourier-plane combination. In this paper, we further discuss and develop data combination 
techniques, and use the most suitable to merge the Hi data for the LMC. 

In § 2 in this paper, we summarize the existing Hi data sets from the ATCA and 
Parkes telescopes and describe the method for linearly combining synthesis and single-dish 
observations. In § 3, we present the combined data set in the form of spatial-spatial map 
integrated over the full velocity range, velocity channel maps, and spatial-velocity maps at 
forty galactic longitude. In § 4, we confirm the Hi giant and supergiant shell candidates 
selected from the ATCA data using the new fully-sampled data, and then compare the size 
distribution of Hi shells with the size distribution of discrete Hi clouds in the LMC. In § 5, 
we give a summary of this paper. 

2. Observation 

The Australia Telescope Compact Array (ATCA) has been used in mosaic mode to 
survey a region 10° x 12° covering the LMC and at an angular resolution of l!0, corresponding 
to a spatial resolution of 15 pc in the LMC. The detailed observations and data reductions 
are described in Kim et al. (1998). The efficient observing techniques and mosaicing methods 
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with the ATCA are described in Staveley-Smith et al. (1997). 

The observing band was centered on 1.419 GHz, corresponding to a central heliocentric 
velocity coverage from —33 to +627 km s -1 , at a velocity resolution of 1.65 km s" 1 . Imag- 
ing of the observed visibility data has been performed using the standard linear mosaicing 
scheme (Sault, Staveley-Smith, and Brouw 1996). After deconvolution the final cube was 
constructed by convolving the maximum entropy model with a circular Gaussian beam of 
FWHM 1' x 1' and adding the deconvolution residuals back in. The RMS noise of the final 
map, determined from the line-free parts of the final data cube, is 31 mJy beam^ 1 . The 
corresponding brightness temperature sensitivity is 2.5 K, which is equivalent to a column 
density sensitivity of 1.6 x 10 19 cm -2 per 1.649 km s -1 channel. 



2.1. Parkes H I Observations 

To provide the complete map at low spatial resolution to complement the ATCA data, 
observations were taken with the inner 7 beams of the Parkes Multibeam receiver (Staveley- 
Smith et al. 1996) on 1998 December 13 to 17 (Staveley-Smith et al. 2001). The receiver 
was scanned across the LMC in orthogonal east-west and north-south great circles and the 
receiver was continuously rotated such that the rotation angle was always at 19?1 to the scan 
trajectory, thus ensuring uniform spatial sampling of the sky. The area covered was 13° by 
14° in RA and Dec, respectively, and centered on 05 h 20 m , -68°44' (J2000). In a single scan, 
the spacing between adjacent tracks is 9'5, which is smaller than the mean FWHP beam 
width of lA'.l, but greater than the Nyquist interval (X/2D) of 5^7. Therefore, six scans are 
interleaved in each of the principal scan directions, resulting in a final track spacing of V.6. 
In total, 12 x 6 RA scans and 11x6 Dec. scans were made. Seven scans were dropped or 
edited out due to drive problems, leaving a total of 131 scans consisting of a total of 29 hours 
of on-source integration on each of 7 beams. The average integration time per beam area is 
360 seconds for both polarizations. 

The scan rate of the telescope was 1?0 min -1 and the correlator was read every 5 seconds. 
Therefore, the beam was slightly broadened in the scan direction to 14.'5. After averaging 
orthogonal scans, the effective beam width reduces to 14'3. The central observing frequency 
was switched between 1417.5 and 1421.5 MHz, again every 5 seconds. This allowed the 
bandpass shape to be calibrated without spending any time off-source. A bandwidth of 8 
MHz was used with 2048 spectral channels in each of two orthogonal linear polarizations. 
Hi emission from the LMC appeared within the band, at both frequency settings. After 
bandpass calibration, the data from both settings were shifted to a common solar barycentric 
reference frame. The velocity spacing of the multibeam data is 0.82 km s _1 , but the final 
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cube was Hanning-smoothed to a resolution of 1.6 km s _1 . The useful velocity range in 
the final cube (i.e. after excluding frequency side lobes of the LMC and the Galaxy, and 
band-edge effects) is —66 to 430 km s _1 . 

Bandpass calibration, velocity shifting and preliminary spectral baseline fitting were 
all done using the aips ++ LiveData task. Subsequently, the data were convolved onto a 
grid of 4' pixels using a Gaussian kernel with a FWHP of 8'.0. This broadens the effective, 
scan-broadened, beam width of the inner 7 beams from 14'3 to about 16'4. Residual spectral 
baselines were removed by fitting polynomials in the image domain (miriad task contsub). 
The multibeam data were calibrated relative to a flux density for PKS B 1934— 638 of 14.9 
Jy at the observing frequency. The brightness temperature conversion factor of 0.80 K Jy _1 
was established by an observation of S9 (T B = 85 K, Williams 1973). On the same scale, we 
measured a brightness temperature for pointing 416 in the SMC (00 h 47 m 52.6 s , -73°02'19'.'8, 
J2000) of T B = 133 K, compared with the 137 K measured by Stanimirovic et al. (1999). 
The 3% difference is probably due to the different characteristics of the feeds used in the two 
observations, and residual uncertainties in absolute bandpass calibration. The rms noise in 
the line-free region of the cube is 27 mK, which is close to the theoretical value. 

2.2. Combining ATCA and Parkes Observations 

Several techniques are available to combine the interferometric and single-dish data. 
The data can be combined during a joint maximum entropy deconvolution operation. Al- 
ternatively the single-dish data can be used as a 'default' image in a maximum entropy 
deconvolution of the interferometer data. Another possibility is to feather together (a lin- 
ear merging process) the single-dish and interferometer images. Stanimirovic et al. (1999) 
have used an approach where the interferometer and single-dish image are added together 
before deconvolution, and then deconvolution performed with a modified point-spread func- 
tion. Tests by Stanimirovic et al. (1999) found that the different techniques produced quite 
similar results. Given this, and as we already had a deconvolved interferometric cube, and 
as the computational requirements to re-perform a deconvolution of these data is large, we 
have used an image feathering approach. 

The technique we have used is a variant of the approach described by Schwarz & Walker 
(1991). As the Parkes and ATCA images give accurate representations of the LMC at short 
spacings and mid to long spacings respectively, a composite image can be formed by filtering 
out the short spacing data from the ATCA image, and then adding the Parkes image. This 
process is most easily visualized in the Fourier domain as in Figure 1 shows the expected 
amplitude as a function of spatial frequency of a point source for our observations. The 
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Fourier transform of the Parkes images were added to the final images with no weighting 
(i.e. 'natural' weight in interferometric nomenclature). The deconvolved ATCA data was 
also Fourier-transformed, but by down-weighting the lower spatial frequencies such that the 
combined weight of the Parkes and ATCA data was the same as the response to a 1/0 
Gaussian. The miriad task immerge was used. 

Before combining, the Parkes image was interpolated onto the same coordinate grid 
as the ATCA mosaicked image. Also the residual primary beam attenuation remaining in 
the ATCA image was applied to the Parkes image (the mosaicing process we have used 
does not perform full primary beam correction when this would result in excessive noise 
amplification). In order to perform the combination of the Parkes and ATCA observations, 
we need to ensure that the flux calibration between the two data types are consistent. Ideally 
we would like to find the ratio of the flux density of an unresolved point source in the field. 
We have estimated this calibration factor by examining data in the Fourier plane between 21 
and 31 meters - data in this annulus is well measured by the Parkes and mosaicked ATCA 
observations. After tapering the ATCA data to the same resolution as the Parkes data, we 
found a scale factor of 1.3 minimized the L\ difference between the interferometer and single- 
dish Fourier components (the real and imaginary parts of the data where treated as distinct 
measurements in the fitting). Note that the scale factor (and indeed the entire feathering 
process) requires a good estimate of the resolution of the Parkes image. The effective beam 
size of 16.9 arcmin was adopted as it gave a scale factor independent of spatial frequencies 
(Stanimirovic et al. 1999). 

The resolution of the combined Hi image of the LMC is I/O which is the same as for 
the ATCA interferometer map. The RMS noise of the combined map, determined from the 
line-free parts of the final data cube, is ~19 mJy beam -1 . This corresponds to a brightness 
temperature sensitivity of ~ 2.4 K. 

3. DATA PRESENTATION 

To assist the comparison, we display in Figure 2 the individual channel maps from the 
combined Hi ATCA map with the Parkes single dish map. This figure can be compared 
with Figure 2 of Kim et al. (1998). The individual channel maps have a velocity resolution 
of 1.649 km s _1 and cover a velocity range of V ffBi =205 km s _1 to 1/^^=334 km s _1 . 
However, the Hi emission is detected mostly in the velocity range of Vh_bl=190-387 km s^ 1 
. The peak H I surface brightness image and the column density image of the LMC are shown 
in Figure 3 and Figure 4. The peak brightness temperature is 136.7 K at RA=05 h 40 m 43 s , 
DEC = -69°48'49.7" (J2000). The peak column density of 7.3±0.3 xlO 21 cm" 2 in the 
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Vh\el=225— 310 km s 1 assuming that the Hi is optically thin at this point. 

There is a remarkable correspondence between the features of this map and the Hi 
emission obtained from the previous ATCA survey (Kim et al. 1998). For example, the 
spiral structure is clearly seen in the individual channel maps (Figure 2b) and the peak H I 
surface brightness map (Figure 3) as well as its integrated map (Figure 4). In contrast to the 
optical image of the LMC (Kim et al. 1999, Figure 5a), both ATCA map and ATCA+Parkes 
combined map show that the H I distribution is uniform and no bar feature corresponding to 
the optical one. The huge Hi hole of diameter about 1.2 kpc resides between the southern 
(de Vaucouleurs & Freeman 1972) and northern spiral arm. 

The large-scale distribution of hydrogen in 30 Dor region and the south of 30 Dor region 
in the coordinate range 05^ 49 m < RA < 05 h 36 m , -73° 00' < DEC < -68° 30' displays 
two large sheets of gas having a relative difference in the line-of-sight velocity about 40 
km s _1 but with relatively small internal velocity dispersions. The general shape of this 
feature seen in the Position- Velocity (P — V) diagram (Figure 5) can be explained as the 
disk (D-) component and a second surface (L-component), discussed by Luks and Rohlfs 
(1992), which is possibly a region affected by the ram-pressure associated with the motion 
of the LMC through the outer halo of our Galaxy (Kim et al. 1998). However, many arclike 
structures seen in the P — V diagram are likely to originate from expanding shells (Kim et 
al. 1999), although it is very difficult to distinguish between the components from expanding 
shells and intrinsic separate disk components (e.g., the L- and D- components) as a result 
of combined action of multiple expanding shells and local random motion. 



4. SIZE SPECTRUM OF H I SHELLS AND CLOUDS 

The structure of the neutral atomic ISM in the LMC shows a complex distribution of 
H I emission, which is chaotic with hundreds of clouds, shells, arcs, rings, and filaments. The 
shell-like structures seen in the ATCA map still dominate the structure of the neutral atomic 
ISM in the LMC revealed by the ATCA+Parkes combined map. The previous visual survey 
of Hi shell candidates chosen from the ATCA Hi data cube (Kim et al. 1999) has been 
re-investigated from the current combined data set. Here we confirm that the Hi supergiant 
shells reported in Kim et al. (1999) are indeed seen in the combined data sets and presented 
in Figure 6 in this paper. Figure 6 is the same as Figure 2 of Kim et al. (1999). 

The physical parameters of the individual supergiant shells are summarized in Table 
1. The mean column density of neutral hydrogen in the LMC is in the order of 2.8 x 10 21 
cm -2 . Assuming a mean column density is distributed uniformly over the disk thickness of 
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A = 360 pc (Kim et al. 1999), we estimate a mean gas particle density of nn ~ 2 cm -3 . We 
may estimate the amount of kinetic energy of the H I gas associated with the expanding H I 
shells, using their derived sizes and expansion velocities. A total predicted kinetic energy of 
the interstellar gas associated with Hi expanding supershells is 3.9 ± 1.3 x 10 53 ergs over 
the mean dynamical age of the expanding shell ~ 6 Myr. This result is remarkably similar 
to the total kinetic energy deposited from the stellar winds ~ 4.3 x 10 53 ergs over 6 Myr, 
referring the results derived from a total ionizing flux of ~ 6.7 x 10 51 photons s~ x for UV 
sources in the LMC (Smith et al. 1987) and using the relationship between the ionizing 
photon flux and the stellar wind mechanical luminosity, L w /N c = 3.2 x 1CT 13 ergs photon -1 
(Wilson 1983). 

The distribution by number of the Hi shells as a function of their radius is presented 
in Figure 17 of Kim et al. (1999). In the range 100 - 1000 pc, the data are consistent with 
a power-law distribution of slope s = —1.5 ± 0.4. In order to compare the size spectrum 
of the Hi shells in the LMC with a size spectrum of Hi clouds, we investigate the Hi 
cloud candidates from the ATCA+Parkes combined data cube. It is well known that the 
interstellar medium, atomic as well as molecular, is distributed in a hierarchical ensemble of 
clouds. Such clouds or clumps could be a condensation formed during the thermally unstable 
cooling. The gas may initially be in either an atomic or a molecular state depending on the 
local physical conditions (density, excitation temperature, etc.). As the volume densities 
increase, the gas phase turns to molecular. Recent CO(J=l— 0) and CO(J=2— 1) study of 
the two nearly face-on galaxies NGC 628 and NGC 3938 shows that the velocity dispersion is 
remarkably constant with radius, 6 km s _1 for NGC 628 and 8.5 km s -1 for NGC 3938, and 
of the same order as the Hi velocity dispersion (Combes & Becquaert 1997). The similarity 
of the CO and H I dispersions suggests that the two components are well mixed, and are only 
two different phases of the same kinematical gas component. The position of Hi clumps can 
be well matched with CO emission maps (Cohen et al. 1988; Israel et al. 1993; Fukui et 
al. 1999). A large fraction of the Hi clumps have a detectable 100 /im emission. Similarly, 
where the H I clumps are more intense, and the clumps are associated with CO clouds. The 
large-scale association between H I and CO is clearly of central interests for studies of cloud 
and star formation in galaxies (Elmegreen & Elmegreen 1987). 

Shapes of Hi clouds are difficult to define and the identification of clouds is still a 
subjective issue. However, majority of Hi clouds can be redefined as Hi clumps at small 
scales in either larger sheets or filaments abound rather than spherical blobs. The present 
H I aperture synthesis survey of the LMC is particularly well suited to reach the statistical 
characteristics of Hi clouds or clumps. Since the distance of the LMC is known as 55 kpc 
(Feast 1991) and the LMC is nearly face-on disk galaxy so that confusion along the line of 
sight is negligible. We have identified and catalogued Hi clouds in the LMC by defining 
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a cloud to be an object composed of all pixels in right ascension, declination, and velocity 
that are simply connected and that lie above the threshold brightness temperature (Scoville 
et al 1987; Lee et al. 1990; Lee et al. 1997). We applied this method rather than gaussian 
clumping method as the other may generate all gaussian clumps which are not realistic in 
general. In fact, most of clouds do not have gaussian profiles, and especially it is true for HI 
clouds. Ideally, one would like to define clouds with a K threshold temperature. However, 
low threshold temperatures are impractical in view of the noise level in the spectra and 
more importantly because of the blending of adjacent clouds which often occurs in crowded 
regions. We have found H I clouds or clumps using the automatic clump identification code 
with three thresholds of the brightness temperature, T B =16 K (~5xTr M s), 32 K, 64 K. 
Final selection of the Hi clouds has been made with high-temperature thresholds in order 
to reduce the blending of emission from unrelated clouds. 

The distribution by number of the Hi clouds as a function of their size is presented 
in Figure 7. The derived sizes of clouds are distributed in a wide range of scales 20 — 400 
pc. The sizes are computed as the square root of the area. Peak of their size distribution 
of Hi clouds or clumps resides in 20 — 30 pc. However, the effective synthesis beam size 
1' limits the size distribution of the smallest Hi clouds. The data are consistent with a 
power-law distribution with slope of between s=— 0.9 ±0.2 and s—— 1.1 ±0.2. Compared 
to the slope found from the size spectrum of CO clouds in the inner Galaxy, — 2.43±0.05 
(Solomon et al. 1987); Ophiuchus, -2.20±0.06; Rosette, -2.26±0.08; Maddalena-Thaddeus, 
-2.46±0.20; M17, -2.09±0.18 (Elmegreen & Falgarone 1996), the power-law distribution 
of the size spectrum of H I clouds is flatter than the typical size spectrum of CO clouds in 
the Milky Way. 

The measured perimeter P and measured enclosed area A of each Hi cloud in a log- 
log plot (Figure 8) gives a set of points lying along with a slope of D/2 = 0.73 ± 0.1. The 
relation between area and perimeter of each identified H I cloud, P oc A D ^ 2 , can determine the 
fractal dimension .D=1.47 ± 0.2 of a cloud boundary (Vogelaar & Wakker 1994; Williams, 
Blitz, & McKee 2000). The measured fractal dimension of Hi clouds in the LMC is a 
similar dimension, D ±1.4, found in many studies of the molecular ISM (Falgarone et al. 
1991; Williams, Blitz, & McKee 2000). For clouds identified with different thresholds of the 
brightness temperature, the fractal dimension D found from the relation between area and 
perimeter is invariant as shown in Figure 8. 

We find that a previous analysis of the holes and shells in the Hi distribution shows 
the same power-law behavior with that of clumps. The implication of the result is that the 
formation of the filament-like or shell-like structures, as well as clumps could be formed with 
the same power-law for dependence on the mass of gas in the galaxy and the stellar energy 
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feedback. Satisfactory theoretical explanation needs to be developed for this observational 
fact. 



5. Summary 

We present the merged images of the LMC from the H I Parkes multibeam observations 
with the aperture synthesis mosaic made by combining data from 1344 separate pointing 
centers using the Australia Telescope Compact Array (ATCA). The images are constructed 
by a Fourier-plane technique and sensitive to structure on scales of 15 pc (for an LMC 
distance of 55 kpc) upwards. The addition of total power data reveals the plateau of diffuse 
Hi emission. We find that total power data is important in recovering the true source 
structures even linear mosaicing technique recovers the 'missing' short spacings in the uv 
plane (Ekers & Rots 1979; Cornwell 1988). The structure of the neutral atomic ISM in 
the LMC reveals the dumpiness of the Hi distribution over the whole of the LMC. We 
have identified and catalogued Hi clouds in the LMC by defining a cloud to be an object 
composed of all pixels in right ascension, declination, and velocity that are simply connected 
and that lie above the threshold brightness temperature. The power-law distribution of the 
size spectrum of H I clouds is similar to the typical size spectrum of the holes and shells in 
the Hi distribution. 

We are indebted to the other team members of ATCA HI mosaic project Mike Kesteven, 
and Dave McConnell. We thank Bruce Elmegreen and Martin White for interesting discus- 
sions. We appreciate anonymous referee for improvement of this paper. 
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Fig. 1. — A representation of the feathering process in the Fourier domain. This gives the 
expected visibility function of a point source from the Parkes and ATCA observations. The 
feathering process filters out the short spatial frequencies of the ATCA data, so that the sum 
of this and the Parkes data give a good representation of the object at all spatial frequencies. 
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Fig. 2. — The individual channel maps for the Hi datacube in the LMC. The heliocentric 
velocity is marked at the top left in each panel. Each panel is the average of five adjacent 
channels of width 1.649 km s -1 giving a panel spacing of 8.2 km s" 1 . The pots cover most 
of the Hi emission in the LMC from 210—334 km s" 1 . Black represents regions with the 
highest brightness temperatures of 136.7 K; white represents K. 

Fig. 2. — continued. 

Fig. 2. — continued. 

Fig. 2. — continued. 

Fig. 3. — The peak H I surface brightness map for the LMC. The grey scale intensity range is 
to 136.7 K. This map is sensitive to the small Hi clouds with the highest opacity along the 
line of sight. In such regions, the brightness temperature will approach the spin temperature 
of the H I . The image is similar to the ATCA-only image of Kim et al. (1998) and emphasizes 
the filamentary bubbly and flocculent structure of the ISM in the LMC. 

Fig. 4. — The H I column density image of the LMC from the combined ATCA and Parkes 
data. The intensity range is to 7.3±0.3x 10 21 H-atom cm" 2 in Hi column density. 

Fig. 5. — Declination Velocity images of the LMC. Each panel is a slice through the LMC 
data cube at pa=0° at the RA specified at the top of each panel. The RA separation of the 
slices is ~2.0 m or ~ 10 beamwidths, and the width of each slice is 23 s , or ~ 2 beamwidths. 
Only a fraction of the data are shown. The grey scale intensity range is to 136.7 K. 

Fig. 5. — continued. 

Fig. 5. — continued. 

Fig. 5. — continued. 

Fig. 6. — Supergiant shells are overlayed on the peak brightness HI map. 

Fig. 7. — Histograms of H I clump size for three thresholds of the brightness temperature, 
T B =16K, 32K, 64K. 
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Fig. 8. — A log-log plot of the measured perimeter versus the measured area in units of 
pixels of Hi clump for three thresholds of the brightness temperature, Tg=16K, 32K, 64K. 



Table 1: List of positions, radii and heliocentric velocities for the Hi supergiant shells (SGSs) 
identified in the LMC. 
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DEC 
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Velocity 


Velocity 




(J2000) 





(km s- 1 ) 


(km s- 1 ) 


SGS1 


04:58:36 


- 73:33:57 


24.3x22.7 


19.0 


253 


SGS2 


04:58:30 


- 68:39:29 


18.7x18.7 


17.0 


272 
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04:59:41 


- 65:44:43 


26.9x26.9 


15.0 


294 


SGS4 


05:02:51 


- 70:33:15 


36.7x28.3 


23.0 
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- 68:31:55 
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